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WAVELET COMPRESSION 



The present invention relates to a method of 
compressing data. 
5 In particular it concerns a method for speeding up the 

operation of a computer set a task of processing the data 
without sacrificing the fidelity of the data in regions of 
special interest. 

In modern computational mechanics, typically, but not 

10 exclusively using Finite Element techniques, solution 
parameters are tabulated for discrete points in the problem 
domain. A typical large dynamic analysis may have a 
complex 3D geometry modelled with a million* elements, with 
results tabulated over a large number of time steps, thus 

15 creating a large four dimensional (4D) data set. 

Consider for example, the problem of visualising a fan 
blade containment analysis of a large fan gas turbine 
engine. From an enormous 4D data set, a stress engineer 
must select appropriate cross-sections in order to 

20 visualise the most significant features as they change over 
the selected time frame. The choice of such cross-sections 
is usually subjective, based on experience and engineering 
intuition. However in models of very complex components 
which are subject to a range of loading conditions, it may 

25 be impossible to be certain that all the significant 
regions of the stress field have been inspected. 
Typically, the computer used to perform the calculations 
required for the analysis of the data set will be based 
around a high specification server, and only the graphical 

30 information needed for each individual view will be pushed 
down a network to an engineer's local machine. 

The analysis of the data set is performed on a high 
specification multi processor machine, whereas the engineer 
performs post-processing on a local machine. The vast 
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quantities of output analysed data supplied by the high 
specification server to the local machine are difficult to 
manage for the local machine so that in some cases it is 
necessary to view a subset of the output analysed data on 
5 the display of the local machine. In order to present a 
slightly different view of the model, or to present a 
section through the model, requires the high specification 
server to perform further processing in order for the new 
image to be shown on the display of the local machine, 

10 Thus it is difficult to show time-step or through model 
movies in real time. The processing may be performed on 
the remote high specification machine, but this is not an 
ideal situation, because interactive demands have a 
negative impact on the efficiency of the high specification 

15 machine for its other batch jobs. 

However, there are three practical disadvantages to 
this technique: 

• the graphical computation can be time consuming 
wherein each cross-section through the model can take 

20 several minutes to load, 

• selection of views is by trial and error, and 

• interactive demands on the server interfere with its 
performance on the other large finite element 
analyses, which it is processing. 

25 For the purpose of exploring the solution domain, 

faster computer response is desirable. Although this would 
seem to indicate that a computer system with improved 
server capability would provide the solution, in practice 
the better the server, the more complex will be the Finite 

30 Element Analyses submitted to it. That is the number and 
complexity of the problems submitted to the machine for 
analysis will expand to fill the available capacity. The 
result is often no improvement at all. 
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The present invention seeks to provide a real solution 
to this problem by implementing a process the effect of 
which is to spare the server from interactive demands, and 
make better use of the graphics capabilities of the local 
5 machines. A key feature underlying the invention is 
selective use of a lower fidelity model to process the data 
set for regions of the model where reduced resolution is 
acceptable, and the graphical display of a model of 
sufficiently reduced size could be managed by the local 

10 workstation. 

Accordingly the present invention provides a computer 
system programmed to process a large data set including 
means for analysing the data set and means for applying a 
data compression technique to the analysed data set such 

15 that the compressed analysed data set has high fidelity in 
regions of interest and has lower fidelity in regions of 
lesser interest. 

Preferably the data compression technique comprises 
the use of a wavelet compression technique . 

20 Preferably the data compression technique produces 

high fidelity in geometric regions of interest at points in 
time of interest. 

Preferably the geometric region has a rapid change in 
the stress field or has a high deformation rate. 

25 Preferably the analysed data set comprises a 4D data 

set . 

Preferably the analysed data set comprises an analysed 
data set of a fan blade containment analysis of a casing. 

Preferably the means for analysing the data set 
30 comprises a means for finite element analysis. 

Preferably the computer system comprises a local 
workstation and a graphical display is produced at the 
location workstation . 
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The present invention also provides a method for 
processing a large data set comprising analysing the data 
set and applying a data compression technique to the 
analysed data set such that the compressed analysed data 
5 set has high fidelity in regions of interest and has lower 
fidelity in regions of lesser interest. 

Preferably the data compression technique comprises 
using a wavelet compression technique. 

Preferably the data compression technique produces 
10 high fidelity in geometric regions of interest at points in 
time of interest. 

Preferably the geometric region has a rapid change in 
the stress field or has a high deformation rate. 

Preferably the analysed data set comprises a 4D data 

15 set. 

Preferably the analysed data set comprises an analysed 
data set of a fan blade containment analysis of a casing. 

Preferably the analysing of the data set comprises 
finite element analysis. 
20 Preferably the method comprises producing a graphical 

display at a local workstation. 

The present invention also provides a computer system 
programmed to process a large data set including means for 
applying a data compression technique to the data set and 
25 means for analysing the compressed data set such that the 
analysis has high fidelity in regions of interest and has 
lower fidelity in regions of lesser interest. 

The present invention further provides a method for 
processing large data sets comprising applying a data 
30 compression technique to the data set and analysing the 
compressed data set such that the analysis has high 
fidelity in regions of interest and has lower fidelity in 
regions of lesser interest. 
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In computational modelling, for example computational 
mechanics, mathematical expressions, e.g. differential 
equations, are solved using finite element analysis, 
computational fluid dynamics, partial difference, boundary 
5 integral method, or other techniques to produce a solution 
given as a set of points within a problem domain. Thus the 
computer system, high specification server, is programmed 
to perform engineering numerical analysis, e.g. finite 
element analysis (FEA), computational fluid dynamics (CFD) , 

10 partial difference, boundary integral method, etc, to 
generate solutions to mathematically posed problems on a 
set of points in the domain of the problem and the solution 
is a set of points within the problem domain, this is the 
analysed data set. The solution is usually arranged in a 

15 grid or finite element mesh. 

The gridded, or meshed, analysed data set is treated 
like an image, which is normally in 4D for geometric 
features modelled over time, but may be more than 4D. 

The high specification server performs an additional 

20 step as part of the batch processing of the analysis. In 
this final step, the high specification server creates a 
substantially compressed version of the analysed data set, 
which can be rapidly processed by the local machine. 

The gridded, or meshed, analysed data set is then 

25 simplified by data compression by applying wavelet 
compression techniques to the gridded, or meshed, analysed 
data set. Thus the wavelet compression technique provides 
a compressed analysed data set, which has high fidelity in 
the regions of interest and has lower fidelity in the 

30 regions of lesser interest. This compressed analysed data 
set is then sent to one or more local workstations where 
the compressed analysed data set is viewed on a display by 
an engineer. The compressed analysed data set is more 
easily manipulated by the workstation to produce different 
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cross-sections views without using the high capability 
server . 

The most significant cross-sectional views are 
presented automatically to the engineer by selecting the 
5 regions which have the most rapid changes in the stress 
field, the highest stress field, the highest deformation 
rates or other variables of interest by selecting stress, 
deformation rate or other variables above a threshold, or 
by using principal component analysis of the analysis data 

10 set or Kriging. 

Wavelets have been used with great success for 2D (two 
dimensional) image compression, with compressions to 5% of 
original data still giving visually acceptable images. 
Thus it may be expected that for 4D (four dimensional) 

15 compression, a reduced model of 0.25% of the original size 
is practical. The wavelets method of compression has the 
advantage over other methods of image compression in that 
sharp contrasts are preserved. It is therefore expected 
that areas which would be of interest to a stress engineer, 

20 for example rapid changes in the stress field, high 
deformation rates, etc, would be preserved in high 
fidelity. The fidelity of the less interesting regions of 
the model would be sacrificed. For example, in a fan blade 
containment analysis, most of the casing is not in the 

25 direct line of impact, and much of it experiences very 
little stress in comparison with the impact sites. 

The region of interest is usually the region where the 
impact takes place, the rest of the geometry is necessary 
to provide the appropriate compliance of the impact region, 

30 but the stress data generated in the regions away from the 
impact region is usually of little consequence, with the 
particular exceptions of regions at boundaries or 
interfaces, for example mounting lugs and bolts. 
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This feature of wavelets may be further harnessed, to 
give a basis on which some of the most significant cross- 
sectional views may be automatically computed and presented 
to the engineer at the beginning of the post-processing 
5 stage. 

The wavelet transform is extremely well suited to 
describing transient signals. The wavelet transform 

analyses in terms of scale, and employs scaled and 
translated combinations of a so-called ^mother' or 

10 analysing wavelet. The mother (and derived daughter) 
wavelets have compact support, that is they are of limited 
duration or spatial extent, and have zero mean. Wavelets 
exhibit the desirable property of being well localised not 
only in time, but also in frequency (or more precisely 

15 characteristic scale) . May different forms of wavelet may 
be selected. Differing families of wavelets may be 

selected on the basis of the trade-off between a compact 
localisation in space and the smoothness of the function. 
The simplest wavelet is the Haar family, which essentially 

20 consists of scaled square step functions. The Coif let 
family resemble derivative-of -Gaussian functions. The 
commonly-used Daubechies family possess fractal 

characteristics. All are governed by strict constraints, 
which guarantee sensible mathematical behaviour. 

25 The choice of wavelet type is likely to influence the 

ratio between the fidelity levels in the region of interest 
of the data set and the fidelity levels in the less 
interesting region of the data set. The choice of wavelet 
enables the high fidelity region of interest to be kept 

30 above a controlled threshold and to minimise the data set 
required or to have a more faithful representation of the 
less interesting regions. 

In an example data was developed from an analytical 
expression sampled over 128 x 128 data points. The 
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analytical expression is summation of three Gaussian 
distributions in 2D: 

f=l 

where a, = 0.1, = 0.3, = 0.8, 6, = 0.4, = 0.6, = 

5 0.7, Xj/^ = 50, if/^ = 1000, t//^ = 10, ^, = 200, = 30, (fij = 

200. The domain was defined as 0< x<l and 0< y <1. In this 
example levels of compression of 90%, 95%, 99%, 99.5% and 
99.8% were achieved using 2D wavelet compression, using the 
Coif let family of wavelets. 

10 A reduction to 0.25% is also appropriate for analysis, 

and for post-processing, so that analysis times may be 
expected to be reduced to 0.000625% of the uncompressed 
original processing time. Corresponding reductions in the 
processing times of much bigger analyses may be 

15 contemplated. This is not likely unless there is a lot of 
degeneracy in the problem. 

In an alternative technique wavelets are used to 
transform mathematical expressions, e.g. differential 
equations, such that the variables are then variables in 

20 wavelet space rather than variables in original space, i.e. 
physical variables in physical space e.g. stress, 
deformation rates etc. (Wavelet families are constructed 
such that they can be used for Hilbert transformations) . 

Computational algebra techniques are used to transform 

25 the problem to wavelet space, and the equations would be 
written in the form of one or more series of mathematical 
expressions involving wavelet functions summed to infinity 
(at this point the equations would be exact, with no 
compression) . 

30 The problem thus constructed in wavelet space is then 

truncated, so that it is given in terms of a finite number 
of wavelet terms. (This is what wavelet compression means). 
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This truncated expression is then recast in a form 
which is appropriate to finite element analysis techniques 
(or other appropriate techniques) using computational 
algebra, because the equations would be too complex to 
5 manage by hand. 

The finite element solver is then run to obtain the 
results to the problem, in terms of wavelet space 
variables . 

Then the solution in wavelet space is transferred back 

10 to original variables in original space e.g. physical 
variables in physical space. 

In this case, the analysis data set would be similar 
to that described above, and similar methods could be used 
to display the results. 

15 Application of wavelet compression to the problem 

prior to analysis, allows the finite element analysis 
solution to take place in high fidelity on the most 
significant geometric locations at the most significant 
points in time, while simultaneously reducing the fidelity 

20 elsewhere in the model. For example, in a fan blade 
containment analysis, most of the casing is not in the 
direct line of impact, and much of it experiences very 
little stress in comparison with the impact sites. 

A reduction to 0.25% is also appropriate for analysis, 

25 and for post-processing, so that analysis times may be 
expected to be reduced to 0.000625% of the uncompressed 
original processing time. Corresponding reductions in the 
processing times of much bigger analyses may be 
contemplated. This could mean that pre-processing time, ie 

30 man-hours spent constructing the finite element mesh, could 
be reduced. 



